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Abstract 

In a class of extra dimensional models with a warped metric and a single brane the photon can 
be localized on the brane by gravity only. An intriguing feature of these models is the possibility of 
the photon escaping into the extra dimensions. The search for this effect has motivated the present 
round of precision orthopositronium decay experiments. We point out that in this framework 
a photon in plasma should be metastable. We consider the astrophysical consequences of this 
observation, in particular, what it implies for the plasmon decay rate in globular cluster stars and 
for the core-collapse supernova cooling rate. The resulting bounds on the model parameter exceed 
the possible reach of orthopositronium experiments by many orders of magnitude. 
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I. INTRODUCTION 

Theories with extra dimensions have been very popular in the last decade 

In a large class of such theories the extra-dimensional space is "warped", i.e. 
the metric scales exponentially along one of the additional dimensions (see Eq. ([T]) later). 
The scaling arises naturally as a solution of Einstein's equations in the extra-dimensional 
"bulk" filled with a negative cosmological constant. This solution has the same origin as 
the inflationary solution, a(t) ~ exp(y/Xt / M p i) , but is "aligned" to scale along one of the 
spatial directions, rather than the time direction. This "inflation along a spatial direction" 
can be arranged by introducing one or more domain walls (branes), tuning their tension(s), 
and replacing the positive cosmological constant of inflation with a negative value. In a 
model with two branes, this setup holds promise for solving the hierarchy problem, as 
the exponents can reduce ratios of vastly different scales to relatively modest numbers 
The same reasoning may explain the smallness of the Yukawa couplings 10j. In a model 
with a single brane, this setup offers an alternative to compact ification |8( and, through the 
AdS/CFT (holographic [ll|) connection [jjj], could in fact describe a four-dimensional world 
with a new conformal sector. 

Of course, in any such theory, one faces the problem of explaining why we see only 
four space-time dimensions. One possible line of argument is that the Standard Model 
fields could be dynamically confined to the four- dimensional Minkowsky defect (brane), as 
discussed already two decades ago (RSO) [l|. In models with warping, it is possible that 
the fields are localized to the brane by the metric itself, i.e., by gravity. The localization 
of the graviton by this mechanism in the model with a single warped extra dimension was 
discussed already in a seminal paper (RSII). A scalar field can be similarly localized 13]. 



Gauge fields are not localized in the minimal setup of RSII, but can be localized if the model 



is extended with additional compact extra dimensions 

The states localized in this way act most of the time as "normal" four-dimensional mass- 
less particles. Under some circumstances, however, they can tunnel into the extra dimen- 
sions, "disappearing" from our world. The tunneling can happen if the state is given a 
nonzero mass, or if it is produced as a virtual state with time-like momentum. Not only is 



1 For an expanded set of references, see, e.g. 
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this suggestion intriguing, but, more importantly, potentially experimentally testable. 

We will focus on the possibility of photon tunneling. This effect would manifest itself as 
unexplained missing energy events at an e + e~ collider. The measurements of the Z boson 
resonance provide considerable bounds on the allowed curvature of the extra dimension (see 



later). Another interesting experimental direction that is 



Deing actively pursued is the search 



for an invisible mode in the orthopositronium decay [151 ] . The orthopositronium serves as 
an e + e~ collider with a hermetic detector. Compared to the Z resonance measurement, one 
obviously loses on the center-of-mass energy, but gains considerably on the sensitivity to the 



branching ratio into the invisible mode. The recently published results [15[ find the bound 
Br(oPs — > extra dim) < 4.2 x 10 -7 , with Br(oPs — » extra dim) < 1CT 8 — 1CT 9 expected in 
the future [l^ ]. 

In this paper, we point out that in the same framework photons in plasma (plasmons) 
should also be subject to the invisible decay. Indeed, plasma modifies the photon dispersion 
relation, in a sense providing it with a mass, thereby opening up the decay channel. In what 
follows, we consider the effects of the additional cooling on the cores of low-mass red giants, 
horizontal branch stars, and core-collapse supernovae. The bounds we find on the model 
parameter exceed the possible reach of the orthopositronium by many orders of magnitude. 

II. TUNNELING INTO EXTRA DIMENSIONS: OVERVIEW 

As already mentioned, the existence of the photon mode localized on the positive tension 
brane in the scenarios with warped extra dimension(s) is well established. Following jjjj], 
let us consider a space with the metric 

ds 2 = a(z) 2 {r\ lLV dx p 'dx u — 5ijd6id6j) — dz 2 . (1) 

Here z labels the infinite warped extra dimension, a(z) = exp (— k\z\). At z = we have a 
domain wall (brane) with positive tension. The variables ^ 6 [0, 27ri?j] label n > 1 additional 
compact dimensions, with radii Ri. The fields are assumed to be independent of 0j. 
The action for the electromagnetic field Ac(x, z) in this space is 

VW\d 4 xdzH^-C, (2) 

i=i n 1 

C = -\F CD F CD . (3) 



In Eq. (|2j), Fcd = dcAo — do Ac and A is a constant with mass dimension 1, which will 
be determined later from the requirement that the standard four-dimensional coupling is 
reproduced. The Latin indices are assumed to run over all coordinates, including the extra 
dimensions. 

The equations of motion in vacuum is dc{y/gF CD ) = 0. In the A 5 = gauge this reads 

Q 

V X "dxF^ = d z (a(z) n+2 d z A,)/a(z) n , (4) 
dnivTdpAv) = 0. (5) 

The Greek indices run through 0, 1, 2, 3 (our space-time) and 7]^ is the usual Minkowsky 
metric. As pointed out in [5], this system of equations has an obvious solution that is 
independent of z, A^(x,z) — > A^x). Eq. (jSJ) is trivially satisfied in this case, while Eq. PJ 
becomes the usual Maxwell's equation for a massless photon. This solution describes the 
zero mode localized on the brane. The reason this is so is because the eigenfunctions are 
normalized with the integration measure J dza n (hence the need to introduce the compact 
dimensions) . 

In general, the eigenfunctions are plane waves, e~ tpx , as a function of x = 0, 1, 2, 3, owing 
to the fact that the Poincare invariance along the brane is preserved. The z dependence 
of the eigenfunctions is given by the eigenmodes of the operator on the right hand side of 
Eq. (j3J). Denoting the eigenvalue by m 2 , we find that m 2 = p 2 and 

- d 2 z A,(z) + (2 + n) fcsign(^) d z A,(z) = e 2k ^m 2 A,(z). (6) 

From the four- dimensional point of view, the higher modes behave as massive photons 
(Eq. (Tj0) for a given value of m 2 ^ takes the form of the Proca equation). As we will see 
shortly, their eigenfunctions are strongly suppressed on the brane. 

Further physical insight can be gained by recasting this equation in the Schrodinger form. 
After changing the variable z — > s = sign(z) [exp (k\z\) — 1] and the redefinition of the fields 
as A^(z) — > <p^(z) = A^(z) exp [-k\z\(n + l)/2] we get 



ld?_ (n + l)(n + 3) _ n + 1 
2ds 2 8(|s| + l) 2 2 [S) 



■2 



(7) 



This transformation is similar to what was done in |8j for the graviton. Not only has the 
first derivative term disappeared, but also the measure with which <f) is normalized is trivial, 



f ds. Eq. ([7]) thus describes a non-relativistic Schrodinger problem and the usual physical 
intuition fully applies here. We have a particle of unit mass in a "volcano" potential 2 , with 
a confining 5-function at the origin and a positive barrier outside that slopes off to zero as 
\s\ — > oo. This potential can support a single bound state of zero energy with the wave 
function 0o( s ) = a/^/2(1 + |s|)~^ n+1 ^ 2 , which corresponds to the flat solution of the original 
equation. 

It is clear in this picture that the spectrum of states residing away from the origin starts 
from zero energy and is continuous. This means that the localized state is only marginally 
bound: an infinitely small perturbation to this setup that lifts the zero mode, — > E' = 
m 2 /2k 2 (for example by decreasing in absolute value the coefficient of the 5-function) of 
the localized state makes it metastable. The particle can then tunnel through the potential 
barrier and escape from the brane [5|. The eigenvalue in this case becomes complex and the 
eigenfunction at \z\ — > oo has an asymptotic form of outgoing plane waves. 

The decay rate due to tunneling for this class of problems can be estimated as follows. 

The turning points of the tunneling on either side of the brane are given by the condition 

{n + l)(n + 3)/8 ( | so | + l) 2 = E' . For s > |so| the solution asymptotes to the plane wave, 

a(.E")e _4v/2 ~^ 7 ' s ', while for s < |so| it can be approximated by the unperturbed function, 

A/n/2(l + |s|) _( ' n+1 ^ 2 . The amplitude of the plane wave a(E') (the barrier penetration factor) 

can then be estimated as roughly the unperturbed solution at the turning point, |a(-E')| 2 ~ 

(_E')( n+1 )/ 2 . The flux away from the brane computed at large \s\ equals 2\a(E')\ 2 y/2E' ~ 
(£/)(n+2)/2_ The mtw of the decay rate r , tQ the 

energy of the metastable state E' is 
~ (E') n / 2 , true in any system of units. In the normal units in which the energy is m, we 



thus obtain for the decay rate in the rest frame 17] 



I7 C = c n m(m/k) n . (8) 

The numerical coefficient c n can be found by considering the properties of the exact solution 
given by the Hankel functions [l7|. We find 

c n = (nn)/(2 n+1 T[n/2 + l] 2 ), (9) 



2 Notice that the coefficients in Eq. ((7]) are different from those for the graviton, which has a localized 
solution even for n = 0. 
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where T denotes the gamma function. Numerically, c n = (1, 7r/4, 1/3, vr/32, 1/45, ...) for 
n = (1,2, 3, 4, 5,...). 

We can also now easily see that the continuum modes residing in the bulk are suppressed 
on the brane. Indeed, they have to tunnel to the brane from the outside. This suppresses 
the wave functions by the barrier penetration factor ~ (m/kY n+1 ^ 2 , making the model 
phenomenologically viable for energies <C k. 

Another way to describe the escape into the extra dimensions is by inspecting the propa- 
gator between two points on the brane 17| • As shown in 17| , the Fourier transform of this 
propagator for a massless localized scalar in the Randall- Sundrum background (n = 0) is 

'pH[%/k) 



2k 2 



p + 1 



71 



(r(2)) 



P_ 
2k 



(10) 



where denotes the Hankel function of the first kind and p = where p 2 is the 

square of the four dimensional momentum. This approach makes it very clear that time-like 
virtual particles are also subject to tunneling. For p 2 > 0, up to the overall normalization 
factor, this propagator has a standard Breit-Wigner form [p 2 + ipT]^ 1 with the imaginary 
part giving the decay rate, r(7* — > extra dim) = / A)\/p 2 ~{p 2 /k 2 ). 

The time-like virtual photon is formed, e.g., in e + e~ annihilation at colliders. The bound 
from the measurements of the Z width at LEP 16|] is k > mz^mz/Ar™*) 1 /' 1 , where 
T l y v < 2.0 MeV is the limit on the additional invisible decay width of Z. For n = 2 

n 

this yields [16j k > 17 TeV. Clearly, getting a tighter bound on T tnv would improve the 



bound, which is the idea behind looking for this process in orthopositronium decay. The 
invisible width is T(oPs — > extra dim) ~ c n m p s (m p s /k) n a 4 (one power of a comes from 
the photon vertex, and three more from the wavefunction of oPs at the origin), compared 
to the standard three-photon width, T(oPs — > 37) ~ m p s a 6 18]. One gets a bound 



k > m oPs (c„m oPs / Ar 



iny "jl/n^/n 



oPs) 



m ? s (c n /BR^ 



o7 s ) 1/n «" 2/n - Compared to the LEP bound, 
one trades a factor of m z /m oPs ~ 10 5 for a factor of (Tf v / (m z a 2 BR% u h )) l / n . Properly 
keeping track of all coefficients, one finds for n = 2, k > 0.5TeV with the present 

accuracy BR(oPs — > extra dim) < 4.2 x 10~ 7 . If the bounds on BR(oPs — > extra dim) are 
improved to the 10~ 10 level, the orthopositronium bound for n = 2 would surpass that of 
LEP. Note that for larger n LEP has a bigger advantage, in particular kiEp{n 00) > rriz, 
while k pJn — > 00) > m p s . 
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III. PLASMON DECAY TO EXTRA DIMENSION 



Let us now consider the effect of plasma on the zero mode. Strictly speaking, one needs 
to specify how the electrons in plasma are localized to our brane, for example with a domain 
wall in a new scalar field. While the fine details will be model dependent, the essential 
features can be obtained by assuming the localization "by hand" , with a delta function, in 
the spirit of 

The effective photon Lagrangian, Eq. ([3]), gains an additional term, 
— (1/2) AcTl CD Ae>5(z), where U CD = (in in) is the photon self-energy in plasma, or 
truncated forward scattering matrix element [20(. The presence of this term changes the 
equation of motion to kd c {^/g~F CD ) = -Yi CD A c 5(z). 

Let us describe the main properties of EL First of all, since Uc5 = Else = 0, it has a 
block diagonal form. In addition we can see from the above equations that n 55 plays no 
role in our gauge (A 5 = 0). Let us therefore concentrate on the 4— dimensional part IL^. In 
the hypothesis of isotropic plasma the tensor U^ u is diagonalizable. Because of the gauge 
invariance, one of the eigenvectors is directed along the photon 4-momentum q, and has 
zero eigenvalue. The others define the directions eW of the different physical polarizations, 
and have in general non- vanishing eigenvalues ir^K If we assume parity invariance the two 
transverse modes have the same eigenvalue 7r T , whereas the eigenvalue of the longitudinal 
mode, 7r L , is in general different. 

In the basis spanned by the equation of motion is diagonal. Eq. <^ in the presence 
of plasma generalizes to 



1 + (n + l)(n + 3) _ /n+J. _ \ ] ^ = m^,„ 



2ds 2 8(|s| + l) 2 V 2 2kA 

where (f)^ denotes the components of <p along e^. Treating the plasma term as a per- 
turbation, in the lowest order of perturbation theory we can write that the zero mode is 
lifted by the energy 5E' = (<f) \7c^5{s)/{2kA)\<f) ) = 7rW(n/2)/(2JfeA). To reproduce the 
four- dimensional phenomenology, we write A = kn/2, m 2 = tc^\ The rest of the argu- 
ment proceeds analogously to the vacuum case considered earlier. The bound state becomes 
metastable and the corresponding decay rate into extra dimensions in the rest frame is given 
by 

if W = c n V^(V^/k) n . (12) 
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The quantities tt^ are related to the plasma frequency, V 7rW = £Ma> p i. Here £W is 
in general a function of the photon energy and momentum. Fortunately, for transverse 



photons it can be shown 



20 



21] to be always close to one, 1 < ( T < a/3/2. Moreover, the 



contribution of longitudinal photons to stellar cooling rates in all cases of interest to us can 
be neglected. 

Lastly, for the purpose of computing the cooling rate we need the decay rate of a moving 
plasmon, r£? ®. The latter is related to the one given in Eq. (fl2l) by the Lorenz factor, 



ppZ (i) _ -p. 



^/o; = Cn (C W ) n+2 ^plKi/^)Ki/^) n , 



(13) 



where w is the energy of the plasmon. 



IV. IMPLICATIONS 



A. Astrophysical Bounds 



The energy loss rate per unit volume in a star is computed as (decay rate) x (energy 
oss) x (photon number density), i.e., as the integral of uY^ over the phase space {e.g., 



3), 



Q T = I ^-((u Jpl ) 3 g((uj pl /T) 



(14) 



where g(x) = (£a/£ 2 — -0/ ( ex P ~ an d the subscript T reminds that this is the 
contribution from transverse photons only. 

Let us consider, first, the stars on the Red Giant (RG) branch. For RG stars (at the 
helium flash) the internal temperature is about T ~ 10 8 K and the density p ~ 10 6 gcm -3 . In 
these conditions, the main standard cooling mechanism is the plasmon decay into neutrinos 
(see, e.g., 



20 



22i|). The rate for this decay is 



SM 



pi 



48 7r 2 a u 
where is a renormalization constant whose value is 



(15) 



1 for transverse photons and 



between and 1 for longitudinal photons 2l|], CV = 0.96 is the vector-current coupling con- 
stant and Gp = 1.166 x 10~ 5 GeV~ 2 is the Fermi constant. The contribution of longitudinal 
photons to this cooling is always less than 10% [23J] and we will neglect it in what follows. 
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We will also neglect the longitudinal contribution to the non-standard cooling, since this is 
certainly a conservative assumption. 

Stellar models with the cooling rate in Eq. (I15p are in good agreement with observations 
of globular cluster populations |24|. To maintain this agreement, we need to constrain any 



additional energy loss to not exceed about twice the standard neutrino luminosity 23j . From 
dHJ, (USD and (USD, we find 

Qed T E d c n (£ T ) n ~ 4 



<5sm Tsm u 



2.8 x 10 3 °- 8n c 



( C^pi \ 

V lOkeV J 



M w 








n— 4 , 


lTeV\ 







(16) 



where u = (Cy g 4 ) / (1536 ir 2 a) ~ 1.5 x 10 3 , g — 0.65 is the weak coupling constant, and we 
set Z T = 1. If we impose that this does not exceeds about 2 we find 

k > C T 4 /n uj p1 B l l n (17) 

where for RG stars 

1 / My/ N 



In the nonrelativistic limit, the plasma frequency is given by oj p \ = 28.7eV(Y e p) 1 / 2 (l + (1.0 x 
lO" 6 ^) 2 / 3 )- 1 / 4 [20|, where p is in units of gem and Y e is the electron fraction. To be 
conservative, we take for ( T its largest value and for u v \ its value in the center of the star 



just before helium flash, u p \ ~ 17.8keV23(, corresponding to p ~ 10 6 gcm~ 3 . This choice 
leads to B = 6.2 x 10 28 and to the bounds 

f M \ 3 

k>183M w — -) ~ 1.4 x 10 21 TeV, (n = 1) , 



^ P i 

A; > 13 M w ( — J ~ 5 x 10 6 TeV, (n = 2) , 

V J 

( M \ 1//3 

k> 4.5M w — ~60TeV, (n = 3) , (19) 

V Wpi / 

which are many orders of magnitude stronger than the direct laboratory bounds 3 . 

Next, let us consider the effects of the extra cooling on the supernova SN1987A. It is 
known that after the explosion this anomalous energy loss cannot significantly exceed Qmux — 



3 The bound for n = 1 should be interpreted to mean that the model is excluded for all energies for which 
it is a valid effective description, possibly up to the Planck scale. 
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3 x 10 33 erg 1 cm 3 s 1 (e.g., 20j), which corresponds to the energy released in neutrinos. 



The plasma frequenc y in a SN core is approximately given by lu% ~ Aap 2 /3n ~ (lOMeV) 2 , 



where // ~ 200MeV 25[] is the electron chemical potential. Therefore u p \ < T ~ 30MeV, 
in which case one can approximate g(x) ~ 2/x 3 , giving Qt — 2TujT 3 /Tr 2 . We find again a 
bound as in Eq. (fTTj) . but with 

= 2(cy,)^» /^y ( ^_y 

n 2 Q Ma x V 30MeV / VlOMeV/ ' V ' 

which implies 

k > 6 x 10 14 TeV , (n = 1) , (21) 
fc > 7 x 10 4 TeV , (n = 2) , 
A;>27TeV, (n = 3) , 

where, to be conservative, we have used ( T = 1. 

Finally, a similar argument applies to stars on the horizontal branch (HB). Typically 
HB stars have an av erag e temperature in the Helium core of T ~ 0.8 x 10 8 K and density 



p ~ 0.5 x 10 4 gcm 3 20|, |26|]. This implies, u p \ ~ 1.5keV < T. In this case, an anomalous 
energy loss cannot be larger than Quaxj p — 10 ergg -1 s -1 in order to havegood agreement 



between the predicted and observed number ratio of HB and RG stars [20|. In this case 
we find (C T = 1) B = 1.7 x 10 30 T 3 (cj p i/1.5keV) 2 , where T 8 = T/10 8 K. The corresponding 
bounds are 

k > 1.1 x 10 21 TeV, (n = l), (22) 
k > 1.1 x 10 6 TeV, (n = 2), 
k>9 TeV , (n = 3) . 

B. Implications for Orthopositronium Decay 

The bounds we just found can be directly translated into the value of the branching ratio 
(BR) necessary to have an analogous bound from the orthopositronium experiment: 

„ IYoPs — > extra dim) - /ro„p c 

BR = v — ; '- ~ 1.5 x 10 5 c, ' 



T(oPs -> 3 7 ) ' \ k 

1.5 x 10 5 / m oPs 
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B \ u; p i 



where, in the last step, we set ( T = 1 for simplicity. Approximately this means BR < 
2 x io- 24+L75n , from RG, BR < 2 x 10 _15_1 - 48n , from SN87A, BR < 2 x 1(T 25 + 2 - 8 ™, from 
HB stars. Thus the astrophysical bounds on the allowed branching ratio of oPs to extra 
dimensions for n = 2 are some 14 orders of magnitude stronger than the present sensitivity 
of the oPs experiments. Moreover, the bound from supernova cooling is at least 8 orders of 
magnitude more stringent than the present experiments for any value of n. 

V. CONCLUSIONS 

We have seen that the models in which the photon is gravitationally trapped on the brane 
face significant constraints from astrophysical considerations. The exact constraint depends 
on the number n of extra compact dimensions. For n = 2 or 3 the AdS curvature k is 
constrained to be orders of magnitude above the electroweak scale. For n = 1 the bound 
extends all the way to the Planck scale. For n > 4, the astrophysical bounds are weaker than 
those coming from the LEP measurement of the Z. For any n, the astrophysical bounds 
imply the rate of orthopositronium decay into extra dimensions that is at least eight orders 
of magnitude smaller than the present experimental sensitivity. It is this implication for the 
ongoing and planned experiments that provides the main motivation for our work. 

A detailed discussion of the implications for the models is beyond the scope of this paper. 
Briefly, our bounds do not exclude the models, but provide significant constraints on them. 
One way to keep the scales in the model close to the electroweak scale is by having a large 
(n > 4) number of extra dimensions. Another possibility is to arrange for an additional 
binding mechanism for the photon, besides gravity. The binding energy in the latter case 
needs to significantly exceed the plasma frequency in the proto-neutron star inside a core- 
collapse supernova (~ 10 MeV). 
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